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p , Abstract 

<^ This manuscript will be published as Chapter 5 in Wiley's textbook Mathe- 

matical Tools for Physicists, 2nd edition, edited by Michael Grinfeld from the 
University of Strathclyde. 

I — — ' The chapter provides an introduction to the basic concepts of Algebraic 

Topology with an emphasis on motivation from applications in the physical 
sciences. It finishes with a brief review of computational work in algebraic 

V^ topology, including persistent homology. 
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1 Introduction 

Topology is the study of those aspects of shape and structure that do not de- 
pend on precise knowledge of an object's geometry. Accurate measurements are 
central to physics, so physicists like to joke that a topologist is someone who 
cannot tell the difference between a coffee cup and a doughnut. However, the 
qualitative nature of topology and its ties to global analysis mean that many 
results are relevant to physical applications. One of the most notable areas of 
overlap comes from the study of dynamical systems. Some of the earliest work 
in algebraic topology was by Henri Poincare in the 1890s, who pioneered a qual- 
itative approach to the study of celestial mechanics by using topological results 
to prove the existence of periodic orbits [57 . Topology continued to play an 
important role in dynamical systems with significant results pertinent to both 
areas from Smale in the 1960s J5S]. More recently, in the 1990s computer anal- 
ysis of chaotic dynamics was one of the drivers for innovation in computational 
topology [371 SOI SUET]. 

As with any established subject, there are several branches to topology: 
General topology defines the notion of 'closeness' (the neighborhood of a 
point), limits, continuity of functions and so on in the absence of a metric. 
These concepts are absolutely fundamental to modern functional analysis; a 
standard introductory reference is [2]. Algebraic topology derives algebraic 
objects (typically groups) from topological spaces to help determine when two 
spaces are alike. It also allows us to compute quantities such as the number 
of pieces the space has, and the number and type of 'holes'. Differential 
topology builds on the above and on the differential geometry of manifolds to 
study the restrictions on functions that arise as the result of the structure of 
their domain. This chapter is primarily concerned with algebraic topology; it 
covers the elementary tools and concepts from this field. It draws on definitions 
and material from Chapter ?? on group theory and Chapter ?? on differential 
geometry. 

A central question in topology is to decide when two objects are the same 
in some sense. In general topology, two spaces, A and B, are considered to be 
the same if there is a homeomorphism, /, between them: / : A — » B is a 
continuous function with a continuous inverse. This captures an intrinsic type 
of equivalence that allows arbitrary stretching and squeezing of a shape and 
permits changes in the way an object sits in a larger space (its embedding), 
but excludes any cutting or gluing. So for example, a circle (x 1 + y 2 = 1) is 
homcomorphic to the perimeter of a square and to the trefoil knot, but not 
to a line segment, and a sphere with a single point removed is homeomorphic 
to the plane. One of the ultimate goals in topology is to find a set of quan- 
tities (called invariants) that characterize spaces up to homeomorphism. For 
arbitrary topological spaces this is known to be impossible [57] but for closed, 
compact 2-manifolds this problem is solved by finding the Euler characteristic 



(see p. 15) and orientability of the surface [531126) . 

What is the essential difference between a line segment and a circle? Intu- 
itively it is the ability to trace your finger round and round the circle as many 



times as you like without stopping or turning back. Algebraic topology is the 
mathematical machinery that lets us quantify and detect this. The idea behind 
algebraic topology is to map topological spaces into groups (or other algebraic 
structures) in such a way that continuous functions between topological spaces 
map to homomorphisms between their associated groupsjj 

In Sections [2j [3] and |4j this chapter covers the three basic constructions of 
algebraic topology: homotopy, homology and cohomology theories. Each has a 
different method for defining a group from the structures in a topological space, 
and although there are close links between the three, they capture different qual- 
ities of a space. Many of the more advanced topics in algebraic topology involve 
studying functions on a space, so we introduce the fundamental link between 
critical points of a function and the topology of its domain in Section[5]on Morse 
Theory. The computability of invariants, both analytically and numerically, is 
vital to physical applications so the recent literature on computational topology 
is reviewed in Section [6j Finally, we give a brief guide to further reading on 
applications of topology to physics. 

2 Homotopy Theory 

A homotopy equivalence is a weaker form of equivalence between topological 
spaces than a homeomorphism that allows us to collapse a space onto a lower- 



dimensional subset of itself (as we will explain in Section 2.3), but it captures 
many essential aspects of shape and structure. When applied to paths in a space, 
homotopy equivalence allows us to define an algebraic operation on loops, and 
provides our first bridge between topology and groups. 

2.1 Homotopy of paths 

We begin by defining a homotopy between two continuous functions f,g:X~^ 
Y. These maps will be homotopic if their images f{X), g(X), can be continu- 
ously morphed from one to the other within Y, i.e. there is a parametrized set 
of images that starts with one and ends with the other. Formally this defor- 
mation is achieved by defining a continuous function F : X x [0, 1] — > Y with 

F(x,0) = f(x) and F(x,l) = 9(x)- 

For example, consider two maps from the unit circle into the unit sphere, 
/, g : S 1 — > S 2 . We use the angle 9 E [0, 2tt) to parametrize S* 1 and the 
embedding x 2 + y 2 + z 2 = 1 in R 3 to define points in S 2 . Define f(9) — 
(cos 9, s'm9, 0) to be a map from the circle to the equator and g(9) — (0, 0, 1), a 
constant map from the circle to the north pole. A homotopy between / and g 
is given by F(8,t) — (cos(irt/2) cos 0, cos(-7rf/2) sin#,sin(7rf/2)) and illustrated 
in Fig. [Tl Any function that is homotopic to a constant function, as in this 
example, is called null-homotopic. 



1 A homomorphism between two groups, <fi : G — > H is a function that respects the group 
operation. That is, </>(a • b) = 0(a) * </>(&), for a,b £ G where ■ is the group operation in G and 
* is the group operation in H . 




Figure 1: The function / : S 1 — > S 2 that maps the circle onto the equator of 
the sphere is homotopic to the function g : S 1 — >■ S" 2 that maps the circle to the 



north pole, 
grey. 



Three sections of the homotopy F : S 1 x [0, 1] — > S 2 arc shown in 



When the domain is the unit interval, X = [0,1], and Y is an arbitrary 
topological space the functions / and g are referred to as paths in Y. A space 
in which every pair of points may be joined by a path is called path-connected. 
It is often useful to consider homotopies between paths that fix their endpoints, 
yo and 2/1, say, so we have the additional conditions on F that for all t € [0,1], 
F(0,t) = /(0) =.g(0) =2/o, andF(l.i) = /(l) =5(1) = 2/i- If a path starts and 
ends at the same point, y — j/i, it is called a loop. A loop that is homotopic 
to a single point, i.e. a null-homotopic loop like the one in the example above, 
is also said to be contractible or trivial. A path-connected space in which 
every loop is contractible is called simply connected. So the real line and the 
surface of the sphere are simply connected, but the circle and the surface of a 
doughnut (the torus) are not. 



2.2 The fundamental group 

We are now in a position to define our first algebraic object, the fundamental 
group. The first step is to choose a base point j/o in the space Y and consider 
all possible loops in Y that start and end at 2/0 ■ Two loops belong to the 
same equivalence class if they are homotopic: given a loop / : [0, 1] — > Y with 
/(0) = /(l) = 2/0, we write [/] to represent the set of all loops that are homotopic 
to /. The appropriate group operation [/] * [g] on these equivalence classes is a 
concatenation of loops defined by tracing each at twice the speed. Specifically, 
choose / and g to be representatives of their respective equivalence classes and 
define f * g(x) = /(2a;) when x e [0, |] and f * g{x) — g(2x — 1) when x € [\, 1]. 
Since all loops have the same base point, this product is another loop based at 
j/o . We then simply set [/] * [g] — [f * g\. The equivalence class of the product 



Figure 2: a) Two non-homotopic loops on the torus with the same base point. 
b) A loop homotopic to the concatenation the loops depicted in (a) . c) Another 
loop in the same homotopy class. 

[/ * g] is independent of the choice of / and g because we can re-parameterize 
the homotopies in the same way as we concatenated the loops. Note, though, 
that the equivalence class [/ * g] consists of more than just concatenated loops; 
Fig. [2] depicts an example on the torus. 

The set of all homotopy equivalence classes of loops based at j/o with the 
operation * forms a group with the identity element being the class of null- 
homotopic loops [e] where e(x) = yo, and the inverse of a loop defined to be the 
same loop traced backwards: [/] _1 = [h] where h(x) — /(l — x). This group is 
the fundamental group of Y with base point yo'- %i(Y, yo). 

The operation taking a topological space to its fundamental group is an 
example of a functor. This word expresses the property alluded to in the 
introduction that continuous maps between topological spaces transform to ho- 
momorphisms between their associated groups. The functorial nature of the 
fundamental group is manifest in the following fashion. Suppose we have a con- 
tinuous function / : X — > Y with f(xo) = Vo- Then given a loop in X with 
base point xo we can use simple composition of the loop with / to obtain a loop 
in Y with base point yo- Composition also respects the concatenation of loops 
and homotopy equivalences so it induces a homomorphism between the funda- 
mental groups: 7Ti(/) : 7Ti(X,o:o) — > 7Ti(V, yo)- When the function / : X — >• Y is 
a homeomorphism, it follows that the induced map 7Ti(/) is an isomorphism of 
their fundamental groups. 

Some further elementary properties of the fundamental group are: 

• A simply connected space has a trivial fundamental group, i.e. only the 
identity element. 

• If Y is path-connected, the fundamental group is independent of the base 
point, and we write tti(Y). 

• The fundamental group respects productajof path-connected spaces: tti {X ) 

Y)=ir 1 (X)xir 1 (Y). 



2 The (Cartesian or direct) product of two spaces (or two groups) X X Y is defined by 
ordered pairs (x, y) where x £ X and jSV. 



• The wedge product of two path-connected spaces (obtained by gluing the 
spaces together at a single point) gives a free producfjon their fundamen- 
tal groups: m(X V Y) = ni(X) * 7Ti (Y). 

• The van Kampen theorem shows how to compute the fundamental 
group of a space X = U U V when U, V, and U C\ V are open, path- 
connected subspaces of X via a free product with amalgamation: iri(X) = 
"K\(U) * wi(V)/N, where N is a normal subgroup generated by elements 
of the form iuin^vd) -1 an d 7 is a loop in ni(U H V), ijj and iy are 
inclusion-induced maps from tti(U C] V) to ni(U) and 7Ti(V) respectively. 
See [29] for details. 



2.3 Homotopy of spaces 

Now we look at what it means for two spaces X and Y to be homotopy 
equivalent or to have the same homotopy type: there must be continuous 
functions / : X —$■ Y and g : Y — > X such that fg : Y — > Y is homotopic to 
the identity on Y and gf : X — > X is homotopic to the identity on X. We can 
show that the unit circle S 1 and the annulus A have the same homotopy type 
as follows. Let 

S 1 = {(r,0) I r=l,0G [0,2tt)} and A = {(r,0) | 1 < r < 2,6 G [0,2tt)} 

be subsets of the plane parametrized by polar coordinates. Define / : S 1 — > A to 
be the inclusion map /(l, 0) = (1,0) and let g : A — > S 1 map all points with the 
same angle to the corresponding point on the unit circle: g(r,8) — (1,0). Then 
gf : S 1 — ?> S 1 is given by gf(l, 9) = (1, 0) which is exactly the identity map. The 
other composition is fg : A — >• A is fg(r,8) = (1,0). This is homotopic to the 
identity iA = (r, 0) via the homotopy F(r, 0, i) = (1— t(r— 1), 0). This example is 
an illustration of a deformation retraction: a homotopy equivalence between 
a space (e.g. the annulus) to a subset (the circle) that leaves the subset fixed 
throughout. 

Spaces that are homotopy equivalent have isomorphic fundamental groups. 
A space that has the homotopy type of a point is said to be contractible 
and has trivial fundamental group. This is much stronger than being simply- 
connected: for example, the sphere S 2 is simply connected because every loop 
can be shrunk to a point, but it is not a contractible space. 

2.4 Examples 

Real space K m , m > 1, all spheres S n with n > 2, any Hilbert space, and any 
connected tree (c.f. Chapter ??) have trivial fundamental groups. 

The fundamental group of the circle is isomorphic to the integers under 
addition: 7ri(S' 1 ) = Z. This can be seen by noting that the homotopy class of 

3 The free product of two groups G * H is an infinite group that contains both G and H 
as subgroups and whose elements are words of the form gi/ii<?2^2 ■ • ■■ 




Figure 3: The projective plane, RP 2 is modelled by the unit disk with opposite 
points on the boundary identified. The black loop starting at (0, 0) is homotopic 
to its inverse with the equivalence suggested by the grey loops. 



a loop is determined by how many times it wraps around the circle. A formal 
proof of this result is quite involved — see Hatcher [29] for details. Any space 
that is homotopy equivalent to the circle will have the same fundamental group, 
this holds for the annulus, the Mobius band, a cylinder, and the 'punctured 
plane' R 2 \ (0,0). 

The projective plane RP 2 is a non-orientable surface defined by identifying 
antipodal points on the boundary of the unit disk (or equivalently, antipodal 
points on the sphere). It has fundamental group isomorphic to Z2 (the group 
with two elements, the identity and r which is its own inverse r 2 = id). To see 
this, consider a loop that starts at the center of the unit disk (0, 0), goes straight 
up to (0, 1) which is identified with (0, —1) then continues straight back up to 
the origin. This loop is in a distinct homotopy class to the null-homotopic loop 
but it is in the same homotopy class as its inverse (to see this, imagine fixing 
the loop at (0, 0) and rotate it by 180° as illustrated in Fig. [3]). 

The fundamental group of a connected graph with v vertices and e edges 
(c.f. Chapter ??) is a free group with n generator^] where n = e — (v — 1) is 
the number of edges in excess of a spanning tree. This demonstrates that the 
fundamental group need not be Abelian (products do not necessarily commute). 

The torus T = S 1 x S 1 so 7Ti(T) = Z x Z. More generally, an orientable 
genus-g surfacerl (g > 2) has fundamental group isomorphic to a hyperbolic 
translation group with 2g generators. 

If a space has a finite cell structure, then the fundamental group can be com- 
puted as a free group with relations in an algorithmic manner; this is discussed 
in Section [6l 



4 A free group with one generator, a say, is the infinite cyclic group with elements 
. . . , a -1 , 1, a, a 2 , .... A free group with two generators a, b, contains all elements of the form 
a xl b J1 a 12 b J2 . . . for integers ife, ji. A free group with n generators is the natural generalization 
of this. 

5 Starting with a sphere, you can obtain all closed oriented 2-manifolds by attaching some 
number of handles (cylinders) to the sphere. The number of handles is the genus. 



2.5 Covering spaces 

The result about the fundamental group of a genus-g surface comes from an- 
alyzing the relationship between loops on a surface and paths in its universal 
covering space (the hyperbolic plane when g > 2). Covering spaces are use- 
ful in many other contexts (from harmonic analysis to differential topology to 
computer simulation), so we describe them briefly here. They are simply a 
more general formulation of the standard procedure of identifying a real- valued 
periodic function with a function on the circle. 

Given a topological space X, a covering space for X is a pair (C,p), where 
C is another topological space and p : C — > X is a continuous function onto 
X. The covering map p must satisfy the following condition: for every point 
x e X, there is a neighborhood U of x such that p -1 (E/) is a disjoint union of 
open sets each of which is mapped homeomorphically onto U by p. The discrete 
set of points p~ 1 (x) is called the fiber of x. A universal covering space is one 
in which C is simply-connected. The reason for the name comes from the fact 
that a universal covering of X will cover any other connected covering of X. For 
example, the circle is a covering space of itself with C = S 1 = {z 6 C : \z\ = 1} 
and Pk{z) — z k for all non-zero integers fc, while the universal cover of the 
circle is the real line with pu : R — > S 1 given by pu(t) = exp(i2irt). The point 
z = (1, 0) G S 1 is then covered by the fiber t e Z C K. We illustrate a covering 
of the torus in Fig. [4] 

When X and C are suitably nice spaces (connected and locally path con- 
nected), loops in X based at xq lift to paths in C between elements of the fiber 
of xq. So in the example of S , a path in R between two integers i < j maps 
under py to a loop that wraps j — i times around the circle. 

Now consider homeomorphisms of the covering space h : C — > C that respect 
the covering map, p(h(c)) — p(c). The set of all such homeomorphisms forms 
a group under composition called the deck transformation group. When 
(C,p) is a universal covering space for X, it is possible to show that the deck 
transformation group must be isomorphic to the fundamental group of X . This 
gives a technique for determining the fundamental group of a space in some 
situations; see Hatcher [55] for details and examples. 

2.6 Extensions and applications 



As we saw in the examples of Section |2.4| the fundamental group of an n- 
dimensional sphere is trivial for n > 2, so the question naturally arises how 
we might capture the different topological structures of S n . To generalize the 
fundamental group, we examine maps from an n-dimensional unit cube I n into 
the space X where the entire boundary of the cube is mapped to a fixed base 
point in X, i.e. / : J" — » X, with f(dl n ) = x . Elements of the higher 
homotopy groups Tr n (X, xq) are then homotopy-equivalence classes of these 
maps. The group operation is concatenation in the first coordinate just as we 
defined for one-dimensional closed paths above. The main difference in higher 
dimensions is that this operation now commutes. 
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Figure 4: The universal covering space of the torus is the Euclidean plane 
projected onto the closed surface by identifying opposite edges of each rectangle 
with parallel orientations. The fibre of the base point on the torus is a lattice of 
points in the cover. The two loops on the torus lift to the vertical and horizontal 
paths shown in the cover. The concatenation of these two loops (see Fig. [2] c) 
lifts to the diagonal path in the cover. The deck transformation group for this 
cover is simply the group of translations that preserve the rectangles which is 
isomorphic to Z x Z = 7Ti(T). 



It is perhaps not too difficult to see that ^(S 12 ) = Z, although the multiple 
wrapping of the sphere by a piece of paper cannot be physically realized in 
K 3 in the same way a piece of string wraps many times around a circle. The 
surprise comes with the result that iTk(S n ) is non-trivial for most (but certainly 
not all) fc > n > 2, and in fact mathematicians have not yet determined all the 
homotopy groups of spheres for arbitrary fc and n 29 . Higher-order homotopy 
groups are a rich and fascinating set of topological invariants that are the subject 
of active research in mathematics. 

One application of homotopy theory arises in the study of topological de- 
fects in condensed matter physics. A classic example is nematic liquid crystals 
which are fluids comprised of molecules with an elongated ellipsoidal shape. The 
order parameter for this system is the (time averaged) direction of the major 
axis of the ellipsoidal molecule: n. For identical and symmetrical molecules, the 
sign and the magnitude of the vector is irrelevant, and so the parameter space 
for n is the surface of the sphere with antipodal points identified — topologi- 
cally RP 2 . The existence of non-contractible loops in M.P 2 is associated with 
the existence of topological line defects in configurations of molecules in the ne- 
matic liquid crystal; see Fig. p5J The fact that 7Ti(IRP 2 ) = Z2 is manifest in the 
fact that two defects of the same type can smoothly cancel one another. The 
second homotopy group ^(MP 2 ) = Z, and this is manifest in the existence of 
point defects ("hedgehogs") where the director field points radially away from 
a central point. See Mermin's original article [35] or Nakahara J3HJ for further 
details. 
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Figure 5: A cross-section through a nematic fluid with a line defect that runs 
perpendicular to the page. Each line-segment represents the averaged direction 
of a single molecule. 

3 Homology 

The fundamental group is a useful invariant but it captures only the one- 
dimensional structure of equivalent loops in a space and cannot distinguish 
between spheres of dimensions greater than two, for example. The higher ho- 
motopy groups do capture this structure but are difficult to compute. The 
homology groups provide a way to describe structure in all relevant dimensions, 
but require a bit more machinery to define. This can seem abstract at first, but 
in fact the methods are quite combinatorial and there has been much recent 
activity devising efficient algorithms to compute homological quantities from 
large data sets (see Section [6]). 

There are a number of different formulations of homology theory that give 
essentially the same results for 'nice' spaces (such as differentiable manifolds). 
The two key ingredients are a discrete cell complex that captures the way a 
space is put together, and a boundary map that describes incidences between 
cells of different dimensions. The algebraic structure comes from defining the 
addition and subtraction of cells. 

The earliest formulation of homology theory is simplicial homology, based 
on triangulations of topological spaces called simplicial complexes. This the- 
ory has some short-comings when dealing with very general topological spaces 
and successive improvements over the past century have culminated in the cur- 
rent form based on singular homology and general cell complexes. Hatcher |29j 
provides an excellent introduction to homology from this modern perspective. 
We focus on simplicial homology here since it is the most concrete and easy to 
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adapt for implementation on a computer. The notation used in this section is 
based on that of Munkres [IT] . 

3.1 Simplicial complexes 

The basic building block is an oriented /c-simplex, <r k , the convex hull of k + 1 
geometrically independent points, {xo,x±, . . . ,#&} C R m , with k < m. For 
example, a 0-simplex is just a point, a 1-simplex is a line segment, a 2-simplex 
a triangle, and a 3-simplex is a tetrahedron. We write a k — (xo, x%, . . . ,Xk) 
to denote a fc-simplex and its vertices. The ordering of the vertices defines an 
orientation of the simplex. This orientation is chosen arbitrarily but is fixed, 
and coincides with the usual notion of orientation of line segments, triangles 
and tetrahedra. Any even permutation of the vertices in a simplex gives another 
simplex with the same orientation, while an odd permutation gives a simplex 
with negative orientation. 

An abstract simplicial complex, C, is a collection of oriented simplices 
with the property that a non-empty intersection of two simplices in C must it- 
self be a simplex in C. If the simplicial complex is finite then it can always be 
embedded in W™ 1 for some m (certain complexes with infinitely many simplices 
can also be embedded in finite-dimensional space). An embedded complex is 
called a geometric realization of C. The subset of M. m occupied by the geometric 
complex is denoted by \C\ and is called a polytope or polyhedron. When a 
topological space X is homeomorphic to a polytope, \C\, it is called a trian- 
gulated space and the simplicial complex C is a triangulation of X. For 
example, a circle is homeomorphic to the boundary of a triangle so the three 
vertices a, b, c and three 1-simplices, (ab), (be), (ca) are a triangulation of the cir- 
cle (see Fig. |6J). All differentiable manifolds have triangulations, but a complete 
characterization of the class of topological spaces that have a triangulation is 
not known. Every topological 2 or 3- manifold has a triangulation, but there 
is a (non-smooth) 4-manifold that cannot have a triangulation (it is related to 
the Lie group Eg |52j). The situation for non-differentiable manifolds in higher 
dimensions remains uncertain. 

3.2 Simplicial homology groups 

We now define the group structures associated with a space X that is triangu- 
lated by a finite simplicial complex C. Although the triangulation of a space is 
not unique, the homology groups for any triangulation of the same space are 
identical (see [H]); this makes simplicial homology well-defined. 

The set of all /c-simplices from C form the basis of a free group called the fc-th 
chain group, Ck(X). The group operation is an additive one; recall that — a k 
is just a k with the opposite orientation, so this defines the inverse elements. In 
general, a fc-chain is the formal sum of a finite number of oriented fc-simplices: 
Cfc = J2i a i a i- The coefficients, at, are elements of another group, called the 
coefficient group that is typically the integers, Z, but can be any Abelian 
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Figure 6: The simplicial complex of a triangle consists of one 2-simplex, three 
1-simplices and three vertices (O-simplices). 



group such as the integers mod 2 Z2, the rational Q, or real numbers R. If the 
coefficient group G needs to be emphasized we write Ck(X; G). 

When k > 1 the boundary operator dk : Ck —> Ck-i maps a fc-simplex 
onto the sum of the (fc — l)-simplices in its boundary. If a k = (x , x%, . . . , x^) 
is a fc-simplex, we have 



dk(cr k ) = Yli-lYixo, ...,Xi,...,X k ) 



i=0 



where (xq, . . . , Xi, . . . , Xk) represents the (fc — l)-simplcx obtained by deleting 
the vertex Xi- The action of the boundary operator on general fc-chains is 
obtained by linear extension from its action on the fc-simplices: dfc(X)j a i a i) = 
Yli a idk{vi)- For fc = the boundary operator is defined to be null: 9o( c o) = 
0. We drop the subscript from the boundary operator when the dimension is 
understood. 

As an example, consider the simplicial complex consisting of a triangle and 
all its edges and vertices, as shown in Fig. |6] The boundary of the 2-simplex 
(a, b, c) is 

d((a,b,c)) = (b,c) - (a,c) + (a,b), 

and the boundary of this 1-chain is: 

d((b, c) - (a, c) + (a, b)) = (c - b) - (c - a) + (b - a) = 0. 

This illustrates the fundamental property of the boundary operator, namely 
that 

d k d k+1 =0 

We now consider two subgroups of Ck that have important geometric in- 
terpretations. The first subgroup consists of fc-chains that map to zero under 
the boundary operator. This group is the group of cycles denoted Zk, it is 
the kernel (or null space) of dk and its elements are called fc-cycles. From the 
definition of do we see that all 0-chains are cycles so Z Q = Cq. The second 
subgroup of Ck is the group of fc-chains that bound a (fc 4- l)-chain. This is the 
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group of boundaries B k , it is the image of dk+i- It follows from (3.2 1 that 
every boundary is a cycle, i.e. the image of dk+i is mapped to zero by dk so 
Bk is a subgroup of Zk ■ In our example of the triangle simplicial complex, we 
find no 2-cycles and a single 1-cycle, (b, c) — (a,c) + (a, b). This 1-cycle is also 
the only 1-boundary so Z\ = B\, and the O-boundaries B are generated by the 
two 0-chains: c — b and c — a (the third boundary a — b — (c — b) — (c — a)). 

Since Bk C Zk, we can form a quotient group, Hk — Zk/ Bk] this is the pre- 
cisely the fe-th homology group. The elements of Hk are equivalence classes 
of fc-cycles that do not bound any k + 1 chain so this is how homology charac- 
terizes /c-dimcnsional holes. Formally, two fc-cycles w, z £ Zk are in the same 
equivalence class if w — z £ Bk] such cycles are said to be homologous. We write 
[z] £ Hk for the equivalence class of cycles homologous to z. For the simple 
example of the triangle simplicial complex, we have already seen that Z\ — B\ 
so that Hi — {0}. The 0-cycles are generated by {a,b,c} and the boundaries 
by {(c— b), (c— a)}, so H has a single equivalence class, [c], and is isomorphic 
toZ. 

The homology groups for some familiar spaces are: 

• Real space W 1 has H = Z and H k = {0} for k > 1. 

• The spheres have H (S n ) = Z, H n (S n ) = Z and H k (S n ) = {0} for all 
other values of k. 

• The torus has H = Z, Hi = Z © Z, H 2 = Z, # fc = {0} for all other fc. 
The 2-cycle that generates H^ is the entire surface. This is in contrast to 
the second homotopy group for the torus, which is trivial. 

• The real projective plane, ~RP 2 has Hq = Z, Hi = Z2, and Hk = {0} for 
k > 2. The fact that H 2 is trivial is a result of the surface being non- 
orientable; even though MP 2 is closed as a manifold, the 2-chain covering 
the surface is not a 2-cycle. 

The combinatorial nature of simplicial homology makes it readily computable. 
We give the classical algorithm and review recent work on fast and efficient 
algorithms for data in Section [6] 

3.3 Basic properties of homology groups 

In general, the homology groups of a finite simplicial complex are finitely gen- 
erated Abelian groups, so the following theorem tells us about their general 
structure (see Theorem 4.3 of |4"I]). 

Theorem 1 If G is a finitely generated Abelian group then it is isomorphic to 
the following direct sum: 

G ~ (Z © • • • © Z) © Z tl © • • • © Z tm . 
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The number of copies of the integer group Z is called the Betti number (3. 
The cyclic groups Z ti are called the torsion subgroups and the U are the torsion 
coefficients; they are defined so that ti > 1 and t\ divides £2 which divides £3 
and so on. The torsion coefficients of the homology group H k (C) measure the 
twistedness of the space in some sense. For example, the real projective plane 
has Hi = Z2, because the 2-chain that represents the whole of the surface has 
a boundary that is twice the generating 1-cycle. The Betti number (3 k of the 
fc-th homology group counts the number non-equivalent non-bounding fc-cycles 
and this can be loosely interpreted as the number of fc-dimensional holes. The 
0-th Betti number, /3 , counts the number of path-connected components of the 
space. 

Some other fundamental properties of the homology groups are as follows: 

• If the simplicial complex has N connected components, X = X 1 U . . . U Xm 
then Ho is isomorphic to the direct sum of N copies of the coefficient group, 
zndH k {X) = H k (X 1 )®---®H k (X N ), 

• Homology is a functor. If / : X — > Y is a continuous function from 
one simplicial complex into another, it induces natural maps on the chain 
groups /j : Ck(X) — > Ck(Y) for each k which commute with the boundary 
operator: <9/j = /j9. This commutativity implies that cycles map to cycles 
and boundaries to boundaries, so that the /j induce homomorphisms on 
the homology groups /* : H k (X) -> H k (Y). 

• If two spaces are homotopy equivalent they have isomorphic homology 
groups (this is shown using the above functorial property). 

• The first homology group is the Abelianization of the fundamental group. 
When X is a path-connected space, the connection between H\(X) and 
■K\{X) is made by noticing that two 1-cycles are equivalent in homology if 
their difference is the boundary of a 2-chain; if we parametrize the 1-cycles 
as loops then this 2-chain forms a region through which one can define a 
homotopy. See [29 for a formal proof. 

• The higher-dimensional homology groups have the comforting property 
that if all simplices in a complex have dimensions < m then H k = {0} 
for k > m. This is in stark contrast to the higher-dimensional homotopy 
groups. 

A particularly pleasing result in homology relates the Betti numbers to an- 
other topological invariant called the Euler characteristic. For a finite sim- 
plicial complex, C, define n k to be the number of simplices of dimension k, 
then the Euler characteristic is defined to be x(C) = ^o — n i + n 2 ■ ■ ■■ The 
Euler- Poincare theorem states that the the alternating sum of Betti numbers 

is the same as the Euler characteristic jJT] : \ = /?o — /3\ + $2 This is one 

of many results that connect the Euler characteristic with other properties of 
manifolds. For example, if M is a compact 2-manifold with a Riemannian met- 
ric, then the Gauss-Bonnet theorem states that 27r% is equal to the integral of 
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Gaussian curvature over the surface plus the integral of geodesic curvature over 
the boundary of M [30]. Further, if M is oricntable and has no boundary then 
it must be homeomorphic to a sphere with g handles and x = 2 — 2g where g is 
the genus of the surface. When M is non-orientable without boundary, then it 
is homeomorphic to a sphere with r cross-caps and x — 2 — r. 

The Euler characteristic is a topological invariant with the property of 
inclusion- exclusion: If a triangulated space X = A U B where A and B are 
both subcomplexes, then 

X (X) = x(A)+x(B)-x(AnB). 

This means the value of % is a localizable one and can be computed by cutting up 
a larger space into smaller chunks. This property makes it a popular topological 
invariant to use in applications |35| . A recent application that exploits the local 
additivity of the Euler characteristic to great effect is target enumeration in 
localized sensor networks [SJ [T7] . The Euler characteristic has also been shown 
to be an important parameter in the physics of porous materials [50] • 

The simple inclusion-exclusion property above does not hold for the Betti 
numbers since they capture global aspects of the topology of a space. Relating 
the homology of two spaces to their union requires more sophisticated algebraic 
machinery that we review below. 

3.4 Homological algebra 

Many results and tools in homology are independent of the details about the 
way the chains and boundary operators are defined for a topological space; they 
depend only on the chain groups and the fact that dd = 0. The study of such 
abstract chain complexes and transformations between them is called homo- 
logical algebra and is one of the original examples in category theory |41j . 

An abstract chain complex is a sequence of Abelian groups and homomor- 
phisms 

• • • -¥ Cfc+i -¥ C k — ^ • • • —^ Co — ► 0, with d k d k+1 = 0. 

The homology of this chain complex is Hk(C) = ker d^/im dk+i- In certain 
cases (such as the simplicial chain complex of a contractible space) we find that 
im dfe + i = ker d k for k > 1 so the homology groups are trivial. Such a sequence 
is said to be exact. 

This property of exactness has many nice consequences. For example, with 
a short exact sequence of groups 

Q^A^B^C^O. 

the exactness means that / is a monomorphism (one-to-one) and g is an epimor- 
phism (onto). In fact g induces an isomorphism of groups C w B/f(A), and 
if these groups are finitely generated Abelian then the Betti numbers satisfy 
j3{B) — /3(C) + (3(A) (where we have replaced f(A) by A since / is one-to-one). 
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Now imagine there is a short exact sequence of chain complexes, i.e. 

-> A k A B k -2i> C k -> 0. 

is exact for all A; and the maps //., g^ commute with the boundary operators 
in each complex (i.e. dsfk = fk-id-A, etc.). Typically, the f k will be induced 
by an inclusion map (and so be monomorphisms), and g k by a quotient map 
(making them cpimorphisms) on some underlying topological spaces. The zig- 
zag lemma shows that these short exact sequences can be joined together into 
a long exact sequence on the homology groups of A, B and C: 

• •• -► H k {A) A H k (B) -^ H k {C) A H k _ x {A) -h H k _ x {B) -► •• • 

The maps /„ and g* are those induced by the chain maps / and g and the bound- 
ary maps 8 are defined directly on the homology classes in H k (C) and H k _i(A) 
by showing that cycles in C k map to cycles in A k _ x via g k , the boundary <9g, 
and fk-i- Details are given in Hatcher [25] , 

One of the most useful applications of this long exact sequence in homology is 
the Mayer- Vietoris sequence, a result that describes the relationship between 
the homology groups of two simplicial complexes, X, Y, their union and their 
intersection. This gives us a way to deduce homology groups of a larger space 
from smaller spaces. 

■ • • A H k (X) ffi H k (Y) -2> H k (X UF)^ 

H k _ x (X n Y) ^4 H k -i[X) © H k ..i(Y) S -^> • • • (1) 
The homomorphisms are defined as follows: 

ife(M) = ([«],-[«]) 

Sk([w], [w 1 ]) = [w + w'] 
v k ([z]) = [dz% 

where in the last equation z is a cycle in X U Y and we can write z — z' + z" 
where z' and z" are chains (not necessarily cycles) in X and Y respectively. 
These homomorphisms are well defined (see, for example, Theorem 33.1 of |41j). 
Exactness implies that the image of each homomorphism is equal to the kernel 
of the following one: im j p = kcr s p , im s p = ker v p , and im v p — ker j p _i. 

The Mayer- Vietoris sequence has an interesting interpretation in terms of 
Betti numbers. First we define N k — ker j k to be the subgroup of H k (X n Y) 
defined by the fc-cycles that bound in both X and in Y. Then by focussing on 
the exact sequence around H k (X U Y) it follows that [13 H] 

/3 k {X U Y) = /3 k (X) + (3 k (Y) ~(3 k (Xr\Y)+ rank N k + rank JV fc _i. 

This is where we see the non-localizable property of homology and the Betti 
numbers most clearly. 
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The Mayer- Vietoris sequence holds in a more general setting than simplicial 
homology. It is an example of a result that can be derived from the Eilenberg- 
Steenrod axioms. Any theory for which these five axioms hold is a type of 
homology theory, see [35] for further details. 

3.5 Other homology theories 

Chain complexes that capture topological information may be defined in a num- 
ber of ways. We have defined simplicial chain complexes above, and will briefly 
describe some other techniques here. 

Cubical homology is directly analogous to simplicial homology using k- 
dimensional cubes as building elements rather than fc-dimensional simpliccs. 
This theory is developed in full in [31j and arose from applications in digital 
image analysis and numerical analysis of dynamical systems. 

Singular homology is built from singular fc-simplices which are simply 
continuous functions from the standard fc-simplex into a topological space X, 
er : (vq, . . . , Vk) — > X. Singular chains and the boundary operator are defined 
as they are in simplicial homology. A greater degree of flexibility is found in 
singular homology since the maps a are allowed to collapse the simplices, e.g. 
the standard A:-simplex, k > may have boundary points mapping to the same 
point in X, or the entire simplex may be mapped to a single point |29j . 

An even more general formulation of cellular homology is made by con- 
sidering general cell (CW) complexes. A cell complex is built incrementally by 
starting with a collection of points X^ C X , then attaching 1-cells via maps of 
the unit interval into X so that end points map into X 1 * ' to form the 1-skeleton 
X^\ This process continues by attaching fc-cells to the (k — l)-skeleton by 
continuous maps of the closed unit fc-ball, <f> : B^ — ¥ X that are homeomorphic 
on the interior and satisfy <f> : dBk — > X^ 1 ' . The definition of the boundary 
operator for a cell complex requires the concept of degree of a map of the k- 
sphere (i.e. the boundary of a (fc + 1) -dimensional ball). For details see Hatcher 

We will see in the section on Morse Theory that it is also possible to define 
a chain complex from the critical points of a smooth function on a manifold. 

4 Cohomology 

The cohomology groups are derived by a simple dualization procedure on the 
chain groups (similar to the construction of dual function spaces in analysis). 
We will again give definitions in the simplicial setting but the concepts carry 
over to other contexts. A cochain <fi k is a function from the simplicial chain 
group into the coefficient group, <p k : Ck(X; G) —> G (recall that G is usually the 
integers, Z, but can be any Abelian group). The space of all fc-cochains forms 
a group called the fc-th cochain group C k (X;G). The simplicial boundary 
operators du ■ Cfc — > C^-i induce coboundary operators 5 k ~ 1 : C k ~ 1 —> C k on 
the cochain groups via 8{<j>) — <j)d. In other words, the cochain 8(<j)) is defined 
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via the action of <j) on the boundary of each fc-simplex a = (xq, Xi, • ■ ■ , x k ): 
5(<f>)(a) = ^2(-iy<j>({x 0> . . . ,Xi, . . . ,Xk)). 

i 

The key property from homology that d^dk+i = also holds true for the 
coboundary: 8 5 = (coboundaries are mapped to zero) so we define the 
fc-th cohomology group as H k {X) = ker <5 fe /im 5 k ~ 1 . Note that cochains 
4> € ker 8 are functions that vanish on the fc-boundaries (not the larger group of 
fc-cycles) , and a coboundary rf G im 8 is one that can be defined via the action 
of some cochain cj) k ^ 1 on the (fc — l)-boundaries. 

The coboundary operator acts in the direction of increasing dimension and 
this can be a more natural action in some situations (such as de Rham cohomol- 
ogy of differential forms discussed below) and also has some interesting algebraic 
consequences (it leads to the definition of the cup product). 

• • ■ ^C k+1 S- C k < <— C° <— 0. 

This action of the coboundary makes cohomology contravariant (induced maps 
act in the opposite direction) where homology is covariant (induced maps act 
in the same direction). If / : X — > Y is a continuous function between two 
topological spaces then the group homomorphism induced on the cohomology 
groups acts as /* : H k (Y) -» H k (X). 

In simplicial homology, the simplices form bases for the chain groups, and 
we can similarly use them as bases for the cochain groups by defining an el- 
ementary cochain a as the function that takes the value one on a and zero 
on all other simplices. For a finite simplicial complex it is possible to represent 
the boundary operator d as a matrix with respect to the bases of oriented k- 
and (k — l)-simplices. If we then use the elementary cochains as bases for the 
cochain groups, the matrix representation for the coboundary operator is just 
the transpose of that for the boundary operator. This shows that for finite 
simplicial complexes, the functional and geometric meanings of duality are the 
same. 

Another type of duality is that between homology and cohomology groups on 
compact closed oriented manifolds (i.e. without boundary). Poincare duality 
states that H (M) = if m _fc(Af) for k € {0, .. . ,m} where in is the dimension 
of the manifold, M; see 29 for further details. 

Despite this close relationship between homology and cohomology on man- 
ifolds, the cohomology groups have a naturally defined product combining two 
cochains and this additional structure can help distinguish between some spaces 
that homology does not. We start with (j> £ C k (X; G) and ip € C l (X; G) where 
the coefficient group should now be a ring R (i.e. R should have both addition 
and multiplication operations; Z, Z p , and Q are rings.) The cup product is 
the cochain <f> ^ if) € C k+l (X;R) defined by its action on a (k + Z)-simplex 
a = (vo, ■ • ■ , Vk+i) as follows: 

(0- ip)(cr) = (t>({v , . . . ,v k ))ip((v k , . . . ,Vk+i)). 
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The relation between this product and the coboundary is: 

8((j> ^ $) = 6<f> ~if> + {-l) k 4> - Sip. 

From this, we see that the product of two cocycles is another cocycle, and if the 
product is between a cocycle and a coboundary, then the result is a coboundary. 
Thus, the cup product is a well defined product on the cohomology groups that is 
anticommutative: [<fi] ^ [ip] — (~l) kl [ip) ^ [<fi] (provided the coefficient ring, G, 
is commutative) . These rules for products of cocycles should look suspiciously 
familiar to those who have read Chapter ??. They are similar to those for 
exterior products of differential forms and this relationship is formalized in the 
next section when we define de Rham cohomology. 

4.1 De Rham cohomology 

One interpretation of cohomology that is of particular interest in physics comes 
from the study of differential forms on smooth manifolds; c.f. Chapter ??. Re- 
call that a differential form of degree fc, to, defines for each point p £ M, an 
alternating multilinear map on fc copies of the tangent space to M at p: 

uj p : T p M x ■ • • x T p M -> R 

The set of all differential k- forms on a manifold M is a vector space, J7 fc (M), and 
the exterior derivative is a linear operator that takes a fc-form to a k + 1-form, 
d k : n k (M) -> n k+1 (M) as defined in Chapter ??. 

The crucial property dd = holds for the exterior derivative. In this context, 
fc-forms in the image of d are called exact, i.e. oj = da for some (fc — 1)- 
form a; and those for which duj = are called closed. We therefore have a 
cochain complex of differential forms and can form quotient groups of closed 
forms modulo the exact forms to obtain the de Rham cohomology groups: 

H k R (M,R) = kerd k /imd k -i 

The cup product in de Rham cohomology is exactly the exterior (or wedge) 
product on differential forms. 

De Rham's theorem states that the above groups are isomorphic to those de- 
rived via simplicial or singular cohomology [5] . And so we see that the topology 
of a manifold has a direct influence on the properties of differential forms that 
have it as their domain. For example, the Poincare Lemma states that if M is a 
contractible open subset of M. n then all smooth closed fc-forms on M are exact 
(the cohomology groups are trivial). In the language of multivariable calculus 
this becomes Helmholtz' theorem that a vector field, V, with curlV = in a 
simply connected open subset of M 3 can be expressed as the gradient of a poten- 
tial function: V = grad/ in the appropriate domain [44]. These considerations 
play a key role in the study of electrodynamics via Maxwell's equations [28 . 
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5 Morse theory 

We now turn to a primary topic in differential topology: to examine the relation- 
ship between the topology of a manifold M and real- valued functions defined on 
M. The basic approach of Morse theory is to use the level cuts of a function 
/ : M — > K and study how these subsets M a = / _1 (— oo, a] change as a is varied. 
For 'nice' functions the level cuts change their topology in a well-defined way 
only at the critical points. This leads to a number of powerful theorems that 
relate the homology of a manifold to the critical points of a function defined on 
it. 

5.1 Basic results 

A Morse function / : M — > R is a smooth real-valued function defined on a 
differentiable manifold M such that each critical point of / is isolated and the 
matrix of second derivatives (the Hessian) is non-degenerate at each critical 
point. An example is illustrated in Fig. [7J The details on how to define these 
derivatives with respect to a coordinate chart on M are given in Chapter ??. 
This may seem like a restrictive class of functions but in fact Morse functions 
are dense in the space of all smooth functions, so any smooth function can be 
smoothly perturbed to obtain a Morse function [32]. Now suppose x £ M is a 
critical point of /, i.e. df(x) — 0. The index of this critical point is the number 
of negative eigenvalues of the Hessian matrix. Intuitively this is the number of 
directions in which / is decreasing: a minimum has index 0, and a maximum 
has index m where m is the dimension of the manifold M. Critical points of 
intermediate index are called saddles since they have some increasing and some 
decreasing directions. 

The two main results about level cuts M a of a Morse function / are that: 

• When [a, b] is an interval for which there are no critical values of / 
(i.e. there is no x £ / _1 ([a, b}) for which df(x) = 0) then M a and M& 
are homotopy equivalent. 

• Let i be a non-degenerate critical point of / with index i, let f(x) = c 
and let e > be such that / _1 [c — e,c + e] is compact and contains no 
other critical points. Then M c+e is homotopy equivalent to M c _ c with an 
i-cell attached. 

(Recall that an i-cell is an i-dimensional unit ball and the attaching map glues 
the whole boundary of the z-cell continuously into the prior space) . The proofs 
of these theorems rely on homotopies defined via the negative gradient flow of 

Gradient flow lines are another key ingredient of Morse theory and allow 
us to define a chain complex and to compute the homology of M. Each point 
x E M has a unique flow line or integral path j x : M — >• M such that 

7,(0) = x and ^^ = -V/( 7 x(*)). 
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Figure 7: Imagine a torus sitting with one point in contact with a plane and 
tilted slightly into the page as depicted. Define a Morse function by mapping 
each point on the torus to its height above the plane. This function has four 
critical points: a minimum, two saddles and a maximum. Left: Five level 
cuts of the height function showing how the topology of a level cut changes 
when passing through a critical point. Right: Gradient flow lines between 
the maximum and the two saddle points, and from each saddle point to the 



minimum. 



Taking the limit as t — > ±00, each flow line converges to a destination and an 
origin critical point. The unstable manifold of a critical point p with index 
i is the set of all x € M that have p as their origin; this set is homeomorphic 
to an open ball of dimension i. Correspondingly, the stable manifold is the 
set of all x that have p as their destination. For suitably 'nice' functions the 
collection of unstable manifolds form a cell complex for the manifold M [4]. 

We can also define a more abstract chain complex which is sometimes referred 
to as the Morse-Smale- Witten complex to reflect the history of its development. 
Let Ci be the chain group derived from formal sums of critical points of index 
i. A boundary operator d : Ci — > Cj_i is then defined by mapping p £ C 8 to a 
sum of critical points ^2 ctjQj £ Cj_i for which there is a flow line with p as its 
origin and q as its destination. The coefficients a.j of the qj in this boundary 
chain are the number of geometrically distinct flow lines that join p and qj (one 
can either count mod 2 or keep track of orientations in a suitable manner). It 
requires some effort to show that 99 = in this setting; see [J] for detailsr] 
Morse homology is the homology computed via this chain complex. 

For finite- dimensional compact manifolds Morse homology is isomorphic to 
singular homology, and we obtain the Morse inequalities relating numbers of 



6 In 2D, think of the flow lines that join a single maximum, minimum pair. In general, such 
a region is bounded by flowlines from the maximum to two saddles and from these saddles to 
the minimum. The boundary of the maximum contains these two saddles and their boundaries 
contain the minimum in oppositely induced orientations. 
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critical points of / : M — > K to the Bctti numbers of M: 

co > Po 

C\ - Co > Pi - Po 
c-i - ci + c > Pi - Pi + p 



J2 {-ir-% = J2 



0<i<m 0<i<m 



where Ci is the number of critical points of / of index i and Pi is the z-th Betti 
number of M. Notice that the final relationship is an equality; the alternating 
sum of numbers of critical points is the same as the Euler characteristic x(M). 
It also follows from the above that c, > Pi for each i. 

5.2 Extensions and applications 

Morse theory is primarily used as a powerful tool to prove results in other set- 
tings. For example, Morse obtained his results in order to prove the existence 
of geodesies on a Riemannian manifold [33] ; most famously, Morse theory forms 
the foundation of a proof due to Smale of the higher-dimensional Poincare con- 
jecture |54j . Morse theory has been extended in many ways that relax conditions 
on the manifold or the function being studied [5] • We mention a few of the main 
generalisations here. 

A Morse-Bott function is one for which the critical points may now not 
be isolated and instead form a critical set that is a closed submanifold. At the 
very simplest level for example, this lets us study the height function of a torus 
sitting flat on a table since the circle of points touching the table is critical [7] . 

The Conley index from dynamical systems is a generalization of Morse 
theory to flows in a more general class than those generated by the gradient 
of a Morse function. For general flows, invariant sets are no longer single fixed 
points but may be periodic cycles or even fractal "strange attractors". In the 
Morse setting, the index is simply the dimension of the unstable manifold of 
the fixed point, but for general flows a more subtle construction is required. 
Conley's insight was that an isolated invariant set can be characterized by the 
flow near the boundary of a neighborhood of the set. The Conley index is then 
(roughly speaking) the homotopy type of such a neighborhood relative to its 
boundary. For details see [HI HS1 [35J • 

Building on Conley's work and the Morse complex of critical points and con- 
necting orbits, Floer created an infinite-dimensional version of Morse homology 
now called Floer homology 4] . This has various formulations which have been 
used to study problems in symplectic geometry (the geometry of Hamiltonian 
dynamical systems) and also the topology of 3- and 4-dimensional manifolds 

M- 

There are a number of approaches adapting Morse theory to a discrete set- 
ting, of increasing importance in geometric modelling, image and data analysis, 
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and quantum field theory. The approach due to Forman is summarized in the 
following section. 

5.3 Forman's discrete Morse theory 

Discrete Morse theory is a combinatorial analogue of Morse Theory for functions 
defined on cell complexes. Discrete Morse functions are not intended to be 
approximations to smooth Morse functions, but the theory developed in [24" ] l2"5 ] 
keeps much of the style and flavor of the standard results from smooth Morse 
theory. 

In keeping with earlier parts of this chapter, we will give definitions for a 
simplicial complex C, but the theory holds for general CW-complexes with little 
modification. First recall that a simplex a is a face of another simplex /3 if 
a C ft in which case j3 is called a coface of a. A function / : C — » M. that 
assigns a real number to each simplex in C is a discrete Morse function if for 
every a^ £ C, f takes a value less than or equal to f{a) on at most one coface 
of a and takes a value greater than or equal to f(a) on at most one face of a. 
In other words, 

#{/3^ +1 > > a | /(ft < f(a)} < 1, 

and 

#{ 7 (p - 1) < a | /( 7 ) > f(a)} < 1, 

where # denotes the number of elements in the set. A simplex a^ p ' is critical 
if all cofaces take strictly greater values and all faces are strictly lower. 

A cell a can fail to be critical in two possible ways. There can exist 7 < a 
such that 7(7) > /(a), or there can exist /3 > a such that f(/3) < f(a). Lemma 
2.5 of [23] shows that these two possibilities are exclusive: they cannot be true 
simultaneously for a given cell a. Thus each non-critical cell a may be paired 
either with a non-critical cell that is a coface of a, or with a non-critical cell 
that is a face of a. 

As noted by Forman (Section 3 of [25 ), it is usually simpler to work with 
pairings of cells with faces than to construct a discrete Morse function on a 
given complex. So we define a discrete vector field V as a collection of pairs 
(a( p \ p( p+1 ') of cells a < ft £ C such that each cell of C is in at most one 
pair of V. A discrete Morse function defines a discrete vector field by pairing 
a^ < /?( p+1 ) whenever f(/3) < /(a). The critical cells are precisely those that 
do not appear in any pair. Discrete vector fields that arise from Morse functions 
are called gradient vector fields. See Fig. [8] for an example. 

It is natural to consider the flow associated with a vector field and in the 
discrete setting the analogy of a flow-line is a ]/-path. A V-path is a sequence 
of cells: 

where (c^, ft) € V, ft > aj+i, and a,; 7^ ttj+i for all i = 0, . . . ,1 — 1. A y-path 
is a non-trival closed V -path if a r — ckq for r > 1. Forman shows that a discrete 
vector field is the gradient vector field of a discrete Morse function if and only 
if there are no non-trivial closed y-paths (Theorem 9.3 of [24]). 
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Figure 8: A simplicial complex with the topology of the torus (opposite edges 
of the rectangle are identified according to the vertex labels) . The arrows show 
how to pair simplices in a gradient vector field. A compatible discrete Morse 
function has a critical 0-cell (a minimum) at a, two critical 1-cells (saddles) at 
edges (b,h) and (d, f) and a critical 2-cell (a maximum) at (e,i,g). 



The four results about Morse functions that we gave earlier all carry over 
into the discrete setting: the homotopy equivalence of level sets away from a 
critical point, adding a critical i-cell is homotopy equivalent to attaching an i- 
cell, the existence of and homology of the Morse chain complex, and the Morse 
inequalities. One of the notable differences between the discrete and continuous 
theories is that flow lines for a smooth Morse function on a manifold are uniquely 
determined at each point, whereas y-paths can merge and split. 

6 Computational topology 

An algorithmic and combinatorial approach to topology has led to significant 
results in low-dimensional topology over the past twenty years. There are two 
main apsects to computational topology: first, research into methods for making 
topological concepts algorithmic, culminating for example, in the beginnings of 
an algorithmic classification of (Haken) 3-manifolds [33] (a result analogous 
to the classification of closed compact 2-manifolds by Euler characteristic and 
orientability). And second, the challenge to find efficient and useful techniques 
for extracting topological invariants from data; see [53] for example. We start 
this section by describing simple algorithms that demonstrate the computability 
of the fundamental group and homology groups of a simplicial complex, and 
then survey some recent advances in building cell complexes and computing 
homology from data. 
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Figure 9: A simplicial complex with the topology of a torus (opposite edges of 
the rectangle are identified according to the vertex labels). A spanning tree T 
with root vertex a is shown in bold. Any closed path that starts and ends at a 
can be decomposed into a sum of loops that lie in T except for a single edge. 

6.1 The fundamental group of a simplicial complex 

In Section [2] we saw that the fundamental group of a topological space could 
be determined from unions and products of smaller spaces or by using a cov- 
ering space. When the space has a triangulation (i.e. it is homeomorphic to a 
polyhedron) there is a more systematic and algorithmic approach to finding the 
fundamental group as the quotient of a free group by a set of relations that we 
summarize below. See [5B] for a complete treatment of this edge-path group. 

Let C be a connected finite simplicial complex. Any path in \C | is nomotopic 
to one that follows only edges in C, and any homotopy between edge-paths 
can be restricted to the 2-simplices of C. This means the fundamental group 
depends only on the topology of the 2-skclcton of C. The algorithm for finding 
a presentation of 7i"i(C) proceeds as follows. 

First find a spanning tree T C C^ i.e. a connected, contractible subgraph 
of the 1-skeleton that contains every vertex of C; see Fig.[9]for an example. One 
algorithm for doing this simply grows from an initial vertex v (the root) by 
adding adjacent (edge, vertex) pairs only if the other vertex is not already in 
T. A non-trivial closed edge-path in C (a loop) must include edges that are not 
in T and in fact every edge in C — T generates a distinct closed path in C^ 1 ' . 
Specifically, for each edge (xi,Xj) 6 C — T there is a closed path starting and 
ending at the root v and lying wholly in T except for the generating edge; we 
label this closed path g.^ . Moreover, any closed path based at v can be written 
as a concatenation of such generating paths where inverses are simply followed 
in the opposite direction: gji = g^- . The gij are therefore generators for a free 
group with coefficients in Z. 

Next we use the 2-skeleton C^ 2 ' to obtain the homotopy equivalences of 
closed edge-paths. Each triangle (xi,Xj,Xk) € C defines a relation in the group 
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via gijgjk9ki = id (the identity) where we also set gij = id if (xi,Xj) £ T. 
Let G(C,T) be the finitely presented group defined by the above generators 
and relations. Then it is possible to show that we get isomorphic groups for 
different choices of T and that G(C,T) is isomorphic to the fundamental group 

M\c\) M- 

If C has many vertices, then the presentation of its fundamental group as 
G(C, T) may not be a very efficient description. It is possible to reduce the num- 
ber of generating edges and relations by using any connected and contractible 
subcomplex that contains all the vertices C^°' C K C C. Generators for the 
edge-path group are labelled by edges in C — K, and the homotopy relations are 
again defined by triangles in C^ 2 \ but we can now ignore all triangles in K. For 
the example of the torus in Fig. [9] we could take K to be the eight triangles in 
the 2x2 lower left corner of the rectangular grid. 

6.2 Smith normal form for homology 

There is also a well defined algorithm for computing the homology groups from 
a simplicial complex C. This algorithm is based on finding the Smith normal 
form (SNF) of a matrix representation of the boundary operator as outlined 
below. 

Recall that the oriented fc-simplices form a basis for the fc-th chain group, Ck- 
This means it is possible to represent the boundary operator, dk '■ Ck — > Ck-i, 
by a (non-square) matrix Af. with entries in {—1,0,1}. The matrix Ak has 
rrik columns and rrik-\ rows where m^ is the number of /c-simplices in C. The 
entry aij is 1 if <x; £ Ck-i is a face of o~j £ Ck with consistent orientation, —1 
if <Xj appears in do~j with opposite orientation and if <7j is not a face of Oj. 
Thus each column of Ak is a boundary chain in Ck-i with respect to a basis of 
simplices. 

The algorithm to reduce an integer matrix to SNF uses row and column 
operations as in standard Gaussian elimination, but at all stages the entries 
must remain integers. The row and column operations correspond to changing 
bases for Ck-i and Ck respectively and the resulting matrix has the form: 



D k = 



B k 




where Bk = 



bx 







Bk is a square matrix with Ik non-zero diagonal entries that satisfy bi > 1 and 
b\ divides 6 2 , divides b 3 , and so on. For a full description of the basic algorithm 
see Munkres [ITj . 

The SNF matrices for dk+i and dk give a complete characterization of the 
fc-th homology group Hk- The rank of the boundary group Bk (im A1-+1) is 
the number of non-zero rows of -D/c+i, i.e., lk+i- The rank of the cycle group 
Zk (ker Ak) is the number of zero columns of Dk, i.e. mk — Ik- The torsion 
coefficients of Hk are the diagonal entries bi of Dk+i that are greater than one. 
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The fcth Betti number is therefore 

(3 k = rank(Z fe ) - rank(_B fe ) = m k - l k - l k+1 . 

Bases for Z k and B k (and hence H k ) are determined by the row and column 
operations used in the SNF reduction but the cycles found in this way typically 
have poor geometric properties. 

There are two practical problems with the algorithm for reducing a matrix 
to SNF as it is described in Munkres [41]. First, the time-cost of the algorithm 
is of a high polynomial degree in the number of simplices; second, the entries 
of the intermediate matrices can become extremely large and create numerical 
problems, even when the initial matrix and final normal form have small integer 
entries. When only the Betti numbers are required, it is possible to do better. 
In fact, if we construct the homology groups over the rationals, rather than 
the integers, then we need only apply Gaussian elimination to diagonalize the 
boundary operator matrices; doing this means we lose all information about the 
torsion however. Devising algorithms that overcome these problems and are fast 
enough to be effective on large complexes is an area of active research. 

6.3 Persistent homology 

The concept of persistent homology arose in the late 1990s from attempts to 
extract meaningful topological information from data [3H1 HZ1 H3] ■ To give a 
finite set of points some interesting topological structure requires the introduc- 
tion of a parameter to define which points are connected. The key lesson learnt 
from examining data was that rather than attempting to choose a single best 
parameter value, it is much more valuable to investigate a range of parameter 
values and describe how the topology changes with this parameter. So persistent 
homology tracks the topological properties of a sequence of nested spaces called 
a filtration ■ • • C C„ C Cj C • • • where a < b 6 I is an index parameter. In a 
continuous setting, the nested spaces might be the level cuts of a Morse function 
on a manifold, so that I is a real interval. In a discrete setting this becomes a 
sequence of subcomplexes indexed by a finite set of integers. In either case as 
the filtration grows, topological features appear and may later disappear. The 
persistent homology group, H k (a, b) measures the topological features from 
C a that are still present in Cf,. Formally, H k (a, b) is the image of the map induced 
on homology by the simple inclusion of C a into Cf,. Algebraically, it is defined 
by considering cycles in C a to be equivalent with respect to the boundaries in 
C b : 

H k {a,b) = Z k (a)/{B k (b)f\Z k {a)). 

Computationally, persistent homology tracks the birth and death of every 
equivalence class of cycle and provides a complete picture of the topological 
structure present at all stages of the filtration. The initial algorithm for doing 
this, due to Edelsbrunner, Letscher and Zomorodian (23) . is surprisingly simple 
and rests on the observation that if we build a cell complex by adding a single 
cell at each step, then (since all its faces must already be present) this cell either 
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creates a new cycle and is nagged as positive, or 'fills in' a cycle that already 
existed and is labelled negative. If a is a negative (fc + l)-cell, its boundary da 
is a fc-cycle and its cells are already nagged as either positive or negative. The 
new cell a is then paired with the most recently added (i.e. youngest) unpaired 
positive cell in da. If there are no unpaired positive cells available, we must 
grow da to successively larger homologous cycles until an unpaired positive cell 
is found. By doing this carefully we can guarantee that a is paired with the 
positive fc-cell that created the homology class of da. Determining whether 
a cell is positive or negative a priori is computationally non-trivial in general 
but there is a more recent version of the persistence pairing algorithm due to 
Zomorodian and Carlsson 13, 53] that avoids doing this as a separate step, and 
also finds a representative fc-cycle for each homology class. 

The result of computing persistent homology from a finite filtration is a 
list of pairs of simplices (a^ k \ r( fe+1 )) that represent the birth and death of each 
homology class in the filtration. The persistence interval for each homology class 
is then given by the indices at which the creator a and destroyer r entered the 
filtration. Some non-trivial homology classes may be present at the final step of 
the filtration, these have an empty partner and are assigned 'infinite' persistence. 
There are a number of ways to represent this persistence information graphically: 
the two most popular techniques are the barcode [TT] and the persistence diagram 
[25] . The barcode has a horizontal axis representing the filtration index; for each 
homology class a solid line spanning the persistence interval is drawn in a stack 
above the axis. The persistence diagram plots the (birth, death) index pair for 
each cycle. These points lie above the diagonal, and points close to the diagonal 
are homology classes that have low persistence. It is possible to show that 
persistence diagrams are stable with respect to small perturbations in the data. 
Specifically, if the filtration is defined by the level cuts of a Morse function on a 
manifold, then a small perturbation to this function will produce a persistence 
diagram that is close to that of the original one [14] . 

6.4 Cell complexes from data 

We now address how to build a cell complex and a filtration for use in persistent 
homology computations. Naturally, the techniques differ depending on the type 
of data being investigated; we discuss some common scenarios below. 

The first construction is based on a general technique from topology called 
the nerve of a cover. Suppose we have a collection of 'good' sets (the sets 
and their intersections should be contractible) U = {U±, . . . , L/jv} whose union 
(J Ui is the space we are interested in. An abstract simplicial complex M(U) 
is defined by making each Ui a vertex and adding a fc-simplex whenever the 
intersection Um D . . . n Uik ^ 0- The nerve lemma states that N(U) has the 
same homotopy type as (J Ui [53] . 

If the data set, X, is not too large, and the points are fairly evenly distributed 
over the object they approximate, it makes sense to choose the Ui to be balls 
of radius a centered on each data p oint : U a = {B(xi, a),Xi € X}. This is often 



called the Cech complex; see Fig. 10 If a < b, we see that N(U a ) C M{Ub) 



29 





Figure 10: Left: Balls of radius o centered on 16 data points. Right: The nerve 
of the cover by balls of radius a gives the Cech complex. In this example the 
complex consists of points, edges, triangles, and a single tetrahedron (shaded 
dark gray). As the radius of the balls increases there are more intersections 
between them and higher-dimensional simplices are created. 



and we have a filtration of simplicial complexes that captures the topology of 
the data as they are inflated from isolated points (a = 0) to filling all of space 
(a —$■ oo). 

A similar construction to the Cech complex that is much simpler to compute 
is the Vietoris-Rips or clique complex. Rather than checking for higher-order 
intersections of balls, we build a 1-skeleton from all pairwise intersections and 
then add a /c-simplex when all its edges are present. This construction is not 
necessarily homotopy equivalent to the union of balls, but is useful when the data 
set comes from a high-dimensional space, perhaps with only an approximate 
metric. 

A drawback of the Cech and Vietoris-Rips complexes is that many unneces- 
sary high-dimensional simplices may be constructed. One way to avoid this is to 
build the Delaunay triangulation. There are many equivalent definitions of this 
widely- used geometric data structure [25]. We start by defining the Voronoi 
partition of space for a data set {x±, . . . , xn} C R m , via the closed cells 

V(xi) = {p such that d(p, Xi) < d(p,Xj) for j ^ i}. 

That is, the Voronoi cell of a data point is the region of space closer to it 
than to any other data point. The boundary faces of Voronoi cells are pieces 
of the (m — l)-dimensional bisecting hyperplanes between pairs of data points. 
The Delaunay complex is the geometric dual to the Voronoi complex: when 
k + 1 Voronoi cells share a (m — fc)-dimensional face there is a fc-simplex in 
the Delaunay complex that spans the corresponding k + 1 data points. See 
Fig. 11 for an example in the plane (m = 2). The geometry of the Voronoi 
partition guarantees that there are no simplices of dimension greater than m in 
the Delaunay complex. 
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Figure 11: Left: The Voronoi diagram of a data set with 16 points. Centre: 
The corresponding Delaunay triangulation. Right: The union of balls of radius 
a centred on the data points and partitioned by the Voronoi cells. The corre- 
sponding triangulation is almost the same as that shown in Fig. 
the tetrahedron there are just two acute triangles. 
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Now consider what happens when we take the intersection of each Voronoi 
cell with a ball centered on the data point, B(xi, a). The Voronoi cells partition 
the union of balls (J B(xi,a) and the geometric dual is a subset of the Delaunay 
complex that is commonly referred to as an alpha complex or alpha shape 
(where alpha refers to the radius of the ball [19, 20 ). By increasing the ball ra- 
dius from zero to some large enough value, we obtain a filtration of the Delaunay 
complex that starts with the finite set of data points and ends with the entire 
convex hull. The topology and geometry of alpha complexes has been used, 
for example, in characterizing the shape of and interactions between proteins 
[2T] . The Betti numbers of alpha shapes are also a useful tool for characterizing 
structural patterns of spatial data [15] such as the distribution of galaxies in the 
cosmic web [55] , 

When the data set is very large, a dramatic reduction in the number of sim- 
plices used to build a complex is achieved by defining landmarks and the witness 
complex. This construction generalises the Voronoi and Delaunay method, so 
that only a subset of data points (the landmarks) are used as vertices for the 
complex, whilst still maintaining topological accuracy. A further advantage is 
that only the distances between data points are required to determine whether 
to include a simplex in the witness complex. See |10) for details, and |12j for an 
extensive review of applications in data analysis. 

Another important class of data is digital images which can be binary 
(voxels are black or white), greyscale (voxels take a range of discrete values), 
or coloured (voxels are assigned a multi-dimensional value). In this setting, the 
structures of interest arise from level cuts of functions defined on a regular grid. 
Morse theory is the natural tool to apply here, although in this application, 
the structures of interest are the level cuts of the function while the domain 
(a rectangular box) is simple. There are a number of different approaches to 
computing homology from such data and this is an area of active research. The 
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works [32 HH1 H] present solutions motivated by applications in the physical 

sciences. 

Guide to further reading 

We give a brief precis of a few standard texts on algebraic topology from math- 
ematical and physical perspectives. 

Allen Hatcher's Algebraic Topology [29] is one of the most widely used texts in 
mathematics courses today and has a strong geometric emphasis. Munkres' [4"T] 
is an older text that remains popular. Spanier [56] is a dense mathematical 
reference and has one of the most complete treatments of the fundamentals 
of algebraic topology. A readable introduction to Morse theory is given by 
Matsumoto [32] and Forman's review article [25] is an excellent introduction to 
his discrete Morse theory. 

Textbooks written for physicists that cover algebraic topology include Naka- 
hara's comprehensive book Geometry, Topology and Physics |43) , Schwarz Topol- 
ogy for Physicists [HI] and Naber Topology, Geometry and Gauge Theory [12] ■ 
Each book goes well beyond algebraic topology to study its interactions with 
differential geometry and functional analysis. A celebrated example of this is 
the Atiyah-Singer index theorem which relates the analytic index of an ellip- 
tic differential operator on a compact manifold to a topological index of that 
manifold, a result that has been useful in the theoretical physics of fundamental 
particles. 

References 

[1] P. Alexandroff and H. Hopf. Topologie. Springer, Berlin, 1935. 

[2] M.A. Armstrong. Basic Topology. Springer- Verlag, 1983. 

[3] C.H. Arns, M.A. Knackstedt, and K.R. Mecke. Reconstructing complex 
materials via effective grain shapes. Physical Review Letters, 91:215506, 
2003. 

[4] Augustin Banyaga and David Hurtubise. Lectures on Morse Homology. 
Kluwer Academic Publishers, Netherlands, 2004. 

[5] Y. Baryshnikov and R. Christ. Target enumeration via Euler characteristic 
integrals. SI AM Journal of Applied Mathematics, 70(9):825-844, 2009. 

[6] P. Bendich, H. Edelsbrunner, and M. Kerber. Computing robustness and 
persistence for images. IEEE Transactions on Visualization and Computer 
Graphics, 16:1251-1260, 2010. 

[7] R. Bott. Nondegenerate critical manifolds. Annals of Mathematics, 
60(2):248-261, 1954. 



32 



[8] R. Bott. Morse theory indomitable. Publications Mathematiques de 
I.H.E.S., 68:99-114, 1988. 

[9] R. Bott and L.W. Tu. Differential forms in algebraic topology. Springer, 
New York, USA, 1982. 

[10] G. Carlsson and V. de Silva. Topological estimation using witness com- 
plexes. In M. Alexa and S. Rusinkiewicz, editors, Eurographics symposium 
on point-based graphics, ETH, Zurich, Switzerland, 2004. 

[11] G. Carlsson, A. Zomorodian, A. Collins, and L. Guibas. Persistence bar- 
codes for shapes. International Journal of Shape Modelling, 11:149-187, 
2005. 

[12] Gunnar Carlsson. Topology and data. Bulletin of the American Mathe- 
matical Society, 46 (2): 255-308, 2009. 

[13] Gunnar Carlsson and Afra Zomorodian. Computing persistent homology. 
Discrete and Computational Geometry, 33:249-274, 2005. 

[14] D. Cohen-Steiner, H. Edelsbrunner, and J. Harer. Stability of persistence 
diagrams. Discrete and Computational Geometry, 37:103-120, 2007. 

[15] Charles C. Conlcy. Isolated invariant sets and the Morse index. AMS, 
Providence, Rhode Island, 1978. 

[16] Charles C. Conley and Robert Easton. Isolated invariant sets and isolating 
blocks. Transactions of the American Mathematical Society, 158:35-61, 
1971. 

[17] J. Curry, R. Ghrist, and M. Robinson. Euler calculus and its applications 
to signals and sensing. Proceedings of Symposia in Applied Mathematics, 
70:75-146, 2012. 

[18] Cecil Jose A. Dclfinado and Herbert Edelsbrunner. An incremental algo- 
rithm for Betti numbers of simplicial complexes. In SCG '93: Proceedings 
of the ninth annual symposium on Computational geometry, pages 232-239, 
New York, 1993. ACM. 

[19] H. Edelsbrunner, D.G. Kirkpatrick, and R. Seidel. On the shape of a set of 
points in the plane. IEEE Transactions on Information Theory, 29(4):551- 
559, 1983. 

[20] H. Edelsbrunner and E.P. Mucke. Three-dimensional alpha shapes. ACM 
Transactions on Graphics, 13:43-72, 1994. 

[21] Herbert Edelsbrunner. Biological applications of computational topology. 
In J.E. Goodman and J. O'Rourke, editors, Handbook of Discrete and Com- 
putational Geometry, chapter 63, pages 1395-1412. CRC Press, Boca Ra- 
ton, FL, 2004. 



33 



[22] Herbert Edelsbrunncr and John Harcr. Computational Topology: An intro- 
duction. American Mathematical Society, Providence, Rhode Island, 2010. 

[23] Herbert Edelsbrunner, David Letscher, and Afra Zomorodian. Topologi- 
cal persistence and simplification. Discrete and Computational Geometry, 
28:511-533, 2002. 

[24] Robin Forman. Morse theory for cell complexes. Advances in Mathematics, 
134:90-145, 1998. 

[25] Robin Forman. A user's guide to discrete Morse theory. Seminaire 
Lotharingien de Combinatoire, 48, 2002. 

[26] G.K. Francis and J.R. Weeks. Conway's ZIP proof. The American Mathe- 
matical Monthly, May 1999. 

[27] P. Frosini and C. Landi. Size theory as a topological tool for computer 
vision. Pattern Recognition and Image Analysis, 9:596-603, 1999. 

[28] P.W. Gross and P.R. Kotiuga. Electromagnetic Theory and Computation: 
A topological approach. Cambridge University Press, Cambridge, UK., 
2004. 

[29] Allen Hatcher. Algebraic Topology. Cambridge University Press, Cam- 
bridge, UK., 2002. 

[30] S.T. Hyde, Z. Blum, T. Landh, S. Lidin, B.W. Ninham, S. Andersson, and 
K. Larsson. The Language of Shape: the role of curvature in condensed 
matter. Elsevier, New York, USA, 1996. 

[31] Tomasz Kaczynski, Konstantin Mischaikow, and Marian Mrozck. Compu- 
tational Homology. Springer- Vcrlag, New York, 2004. 

[32] Yukio Matsumoto. An Introduction to Morse Theory. AMS Bookstore, 
USA, 2002. 

[33] Sergei Matveev. Algorithmic topology and classification of 3-manifolds. 
Springer- Vcrlag, Berlin, 2003. 

[34] Dusa McDuff. Floer theory and low dimensional topology. Bulletin of the 
American Mathematical Society, 43:25-42, 2005. 

[35] Klaus Mccke. Integral geometry in statistical physics. International Journal 
of Modern Physics B, 12(9):861, 1998. 

[36] N.D. Mermin. The topological theory of defects in ordered media. Reviews 
of Modern Physics, 51:591-648, 1979. 

[37] G.B. Mindlin, H.G. Solari, M.A. Naticllo, R. Gilmore, and X.-J. Hou. Topo- 
logical analysis of chaotic time series data from the Bclousov-Zhabotinskii 
reaction. Journal of Nonlinear Science, 1:147-173, 1991. 



34 



[38] K. Mischaikow and M. Mrozek. Conley index. In Bernold Fiedler, edi- 
tor, Handbook of Dynamical Systems, volume 2, chapter 9, pages 393-460. 
Elsevier, New York, NY, USA, 2002. 

[39] M. Morse. The calculus of variations in the large, volume 18 of Colloquium 
Publications. American Mathematical Society, Providence, Rhode Island, 
1934. 

[40] M.R. Muldoon, R.S. MacKay, J.R Hukc, and D.S. Broomhcad. Topology 
from a time series. Physica D, 65:1-16, 1993. 

[41] James R. Munkres. Elements of Algebraic Topology. Addison- Wesley, Read- 
ing, MA, 1984. 

[42] G. Naber. Topology, Geometry and Gauge fields: Foundations. Springer, 
Berlin, 2nd edition, 2011. 

[43] Mikio Nakahara. Geometry, Topology and Physics. Taylor and Francis, 2nd 
edition, 2003. 

[44] C. Nash and S. Sen. Topology and geometry for physicists. Academic Press, 
London, UK, 1983. 

[45] A. Okabe, B. Boots, K. Sugihara, and S.N. Chiu. Spatial Tessellations: 
Concepts and Applications of Voronoi Diagrams. John Wiley & Sons, 
Chichester, 2nd edition, 2000. 

[46] V. Robins. Bctti number signatures of homogeneous Poisson point pro- 
cesses. Physical Review E, 74:061107, 2006. 

[47] V. Robins, J.D. Meiss, and E. Bradley Computing connectedness: An 
exercise in computational topology. Nonlinearity, 11:913-922, 1998. 

[48] Vanessa Robins. Towards computing homology from finite approximations. 
Topology Proceedings, 24:503-532, 1999. 

[49] Vanessa Robins, Peter John Wood, and Adrian P. Sheppard. Theory and 
algorithms for constructing discrete Morse complexes from grayscale digital 
images. IEEE Transactions on Pattern Analysis and Machine Intelligence, 
33(8):1646-1658, 2011. 

[50] C. Scholz, F. Wirncr, J. Gdtz, U. Rile, G.E. Schroder-Turk, K. Mecke, and 
C. Bechingcr. Permeability of porous materials determined from the Euler 
characteristic. Physical Review Letters, 109:264504, 2012. 

[51] A.S. Schwarz. Topology for Physicists. Springer, Berlin, 2002. 

[52] A. Scorpan. The Wild World of ^-Manifolds. American Mathematical 
Society, Providence, Rhode Island, 2005. 



35 



[53] H. Seifert and W. Threlfall. A Textbook of Topology. Academic Press, 1980. 
Published in German 1934. Translated by M.A. Goldman. 

[54] S. Smale. Generalized Poincare's conjecture in dimensions greater than 
four. Annals of Mathematics, 74(2):391-406, 1961. 

[55] S. Smale. Diffcrcntiablc dynamical systems. Bulletin of the American Math- 
ematical Society, 73:747-817, 1967. 

[56] Edwin Spanier. Algebraic Topology. Springer, New York, NY, USA, 1994. 
First published 1966 by McGraw-Hill. 

[57] John Stillwell. Mathematics and Its History. Springer, New York, USA, 
third edition, 2010. 

[58] R. van de Weygaert, G. Vegtcr, H. Edclsbrunner, and et.al. Alpha, Betti 
and the megaparsec universe: on the topology of the cosmic web. Trans. 
Comput. Set, XIV:60-101, 2011. 

[59] Afra Zomorodian. Computational topology. In M. Atallah and M. Blan- 
ton, editors, Algorithms and Theory of Computation Handbook, volume 2 
Special Topics and Techniques, chapter 3. Chapman & Hall/CRC Press, 
Boca Raton, FL, second edition, 2009. 



36 



